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Abstract. Stratospheric Aerosol and Gas Experiment (SAGE 
II) version 7 (v7) ozone profiles are analyzed for their 
decadal-scale responses in the middle and upper stratosphere 
for 1991 and 1992-2005 and compared with those from 
its previous version 6.2 (v6.2). Multiple linear regression 
(MLR) analysis is applied to time series of its ozone num- 
ber density vs. altitude data for a range of latitudes and alti- 
tudes. The MLR models that are fit to the time series data in- 
clude a periodic 1 1 yr term, and it is in-phase with that of the 
1 1 yr, solar UV (Ultraviolet)-flux throughout most of the lat- 
itude/altitude domain of the middle and upper stratosphere. 
Several regions that have a response that is not quite in-phase 
are interpreted as being affected by decadal-scale, dynamical 
forcings. The maximum minus minimum, solar cycle (SC- 
like) responses for the ozone at the low latitudes are similar 
from the two SAGE II data versions and vary from about 5 
to 2.5 % from 35 to 50 km, although they are resolved better 
with v7. SAGE II v7 ozone is also analyzed for 1984-1998, 
in order to mitigate effects of end-point anomalies that bias 
its ozone in 1991 and the analyzed results for 1991-2005 or 
following the Pinatubo eruption. Its SC-like ozone response 
in the upper stratosphere is of the order of 4 % for 1984-1998 
vs. 2.5 to 3 % for 1991-2005. The SAGE II v7 results are also 
recompared with the responses in ozone from the Halogen 
Occultation Experiment (HALOE) that are in terms of mix- 
ing ratio vs. pressure for 1991-2005 and then for late 1992- 
2005 to avoid any effects following Pinatubo. Shapes of their 
respective response profiles agree very well for 1992-2005. 
The associated linear trends of the ozone are not as negative 
in 1992-2005 as in 1984-1998, in accord with a leveling off 
of the effects of reactive chlorine on ozone. It is concluded 
that the SAGE II v7 ozone yields SC-like ozone responses 
and trends that are of better quality than those from v6.2. 


1 Introduction 

Analysis results have been reported for the 1 1 yr solar cy- 
cle (SC) responses (maximum (max) minus minimum (min) 
from 2 to 4 %) and the trends in upper stratospheric ozone 
from long-term satellite measurements (e.g., Dhomse et al., 
2011; McLinden and Fioletov, 2011; Remsberg and Lingen- 
felser, 2010; Fioletov, 2009; McCormack et al., 2007; Ran- 
del and Wu, 2007; Soukharev and Hood, 2006; and Lee and 
Smith, 2003). However, there has been some disagreement 
between the response profiles from stratospheric models and 
from the analyzed data sets (WMO, 2007). Fioletov (2009) 
reported that he had difficulty in obtaining an accurate 1 1 yr 
response for observed upper stratospheric ozone from a time 
series of “merged, partial column ozone data” from the set 
of successive, satellite Solar Backscatter Ultraviolet (SBUV) 
instruments, even though he was able to resolve its 27 day 
ozone response. McLinden and Fioletov (denoted MF, 2011) 
went further and calculated an ozone partial column vs. pres- 
sure response from the Stratospheric Aerosol and Gas Exper- 
iment (or SAGE II) version 6.2 (v6.2) data that agreed with 
that from SBUV, but only after making use of responses and 
trends in temperature that were smaller than indicated by the 
operational temperatures archived along with the SAGE II 
ozone. 

Ozone responses to variations of the UV (Ultraviolet) 
flux ought to be very nearly the same for each solar cycle. 
However, responses diagnosed from the satellite data of the 
past three decades may also contain biases from concurrent 
trends in ozone due to changes in the reactive chlorine, due 
to perturbations from the major volcanic eruptions of 1982 
and 1991, or because of uncertainties for the time series of 
ozone data from the SBUV sensors on successive operational 
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satellites. Randel and Wu (denoted RW, 2007), among others, 
used an equivalent effective stratospheric chlorine (EESC) 
quantity as a proxy for accounting for the ozone trends of 
1979-2005, due to the nonlinear changes in reactive chlorine. 
However, MF cautioned that the effects of EESC changes de- 
pend on the specific units for a given ozone data series, i.e., 
number density vs. altitude or partial column vs. pressure. 
Soukharev and Hood, (denoted SH, 2006) found that simple, 
linear trend terms were adequate for an accounting of the ef- 
fects of the changing chlorine in shorter ozone time series. 
To reduce the perturbing effects from volcanic activity, Lee 
and Smith (2003) excluded data from their ozone series for 
many months following major eruptions. 

Remsberg and Lingenfelser (denoted RL, 2010) diagnosed 
decadal-scale responses from the SAGE II v6.2 ozone time 
series of 1991-2005, and they found good consistency for 
the shape of its response profiles in the upper stratosphere 
vs. those from a representative model and from the Halogen 
Occultation Experiment (or HALOE) ozone. Their SAGE II 
responses were in terms of number density vs. altitude, while 
the ones from HALOE were for mixing ratio vs. pressure. 
However, the response profiles from the model and the ob- 
servations diverged above 50 km. 

The current study is an update of the work of RL, using the 
newly released SAGE II, version 7 (or v7), ozone that is de- 
scribed in Kyrola et al. (2013) and in Damadeo et al. (2013) 
(see also https://eosweb.larc.nasa.gov/project/sage2/sage2_ 
release_v7_notes). Time series analyses are conducted for 
the period 1991-2005, as in RL, and also 1984-1998 to mini- 
mize the effect of end-point anomalies in the ozone following 
the Pinatubo eruption of 1991. In addition, results are shown 
for 1992-2005, in order to assess the effects of the eruption 
on the trends and the analyzed responses. The use of linear 
trends for the three sets of time series is similar to the piece- 
wise linear trend approach used by Kyrola et al. (2013), al- 
though they did not focus on the character of their analyzed 
SC response profiles. 

In the present study proxy terms for the quasi-biennial os- 
cillation (QBO) and the solar cycle are not used, in order 
to be consistent with the previous multiple linear regression 
(MLR) analyses of RL and to show results from just updat- 
ing SAGE II ozone from v6.2 to v7. Section 2 illustrates one 
aspect of the v7 ozone that has led to the improvements in 
the current results. It is also shown that throughout the lati- 
tude/altitude domain, the 1 1 yr or decadal-scale responses are 
almost always in-phase with the solar flux maximum and are 
thereby nearly equivalent to those obtained from a solar UV 
proxy. The SC-like response profiles from v7 and v6.2 are 
compared at low latitudes, along with those from the HALOE 
ozone and from a representative model. Section 3 shows the 
SC-like responses from SAGE II for both 1984-1998 and 
1992-2005, in order to minimize perturbations in the ozone 
following the eruption of Pinatubo. The distributions of as- 
sociated trends in the ozone are also compared for the two 
periods and then related to model simulations of the effects 
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from changes in reactive chlorine that have been reported in 
the literature. The response profiles from SAGE II (number 
density vs. altitude) and from HALOE (mixing ratio vs. pres- 
sure) are compared and related qualitatively for 1992-2005, 
or when the recently reported temperature trends are small 
for the upper stratosphere. Section 4 compares the current 
SAGE II v7 ozone response profiles with previously pub- 
lished results and comments on several of the discrepancies. 
Section 5 is a summary of the findings. 


2 Analysis and results for the decadal-scale responses 
and trends in SAGE II v7 ozone 

Periodic and linear trend terms are applied in the MLR model 
fits of the SAGE II ozone time series for each altitude and lat- 
itude bin. The analyses herein are focused on the max minus 
min responses (in percent) for a periodic 1 1 yr term as a func- 
tion of altitude and latitude for the v7 ozone of 1991-2005, as 
in RL. Responses are obtained at altitudes from 27.5 to 55 km 
at increments of 2.5 km and at twelve latitude bins spaced by 
10° and centered from 55° S to 55° N. The latitude incre- 
ments for each of the bin-averaged, SR and SS points in an 
individual time series have a width of 1 0° for the higher lati- 
tudes (35, 45, and 55°), but widen to 20° at the lower (5, 15, 
and 25°) latitudes to accommodate the fewer samples there. 
A total of 144 ozone time series are analyzed. 

A low altitude limit of 27.5 km was chosen for the anal- 
yses to avoid significant interfering effects of the volcanic 
aerosols on the retrieved ozone and because the period of the 
QBO forcing for ozone is somewhat variable for the lower 
stratosphere. As in RL, proxy terms are not used for the MLR 
models. Instead, Fourier analysis of the time series of the 
residuals after removing the seasonal terms almost always in- 
dicates that there are two significant, interannual terms hav- 
ing periods of order of 28 (QBO-like) and 21 months (sub- 
biennial term denoted as IA). The I A term is the result of the 
difference interaction between the QBO and annual cycles 
(Dunkerton, 2001; Tung and Yang, 1994). Its period varies 
from 20 to 22 months, most likely because of slight varia- 
tions in the period of the actual QBO forcing. Remsberg et 
al. (2001, Fig. 4d) contains an example plot of the diagnostic 
that is used to check for the presence of periodic structure 
in the residual time series, and Remsberg (2008, Figs. 7, 8) 
gives the latitude/pressure distributions of the amplitudes of 
the QBO and IA terms in the HALOE ozone data. Those two 
interannual terms are also present in the SAGE II v7 ozone 
time series, so they are included in the current MLR models 
along with the annual and semiannual terms. Amplitude of 
the I A term is often near to that of the QBO in the subtropi- 
cal middle stratosphere, and its inclusion in the MLR models 
leads to improved fits for the 1 1 yr and trend terms at those 
latitudes. 

Salby et al. (1997) and Baldwin and Dunkerton (1998) 
noted that there is a decadal-scale interaction between the 
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QBO and biennial oscillations that can mimic the solar cy- 
cle, while Smith and Matthes (2008) modeled the response 
of ozone to the SC in the presence of the QBO cycle. How- 
ever, the present analysis simply includes an 1 1 yr periodic 
term as part of the MLR fit to the data. Whenever its maxi- 
mum occurs within one year of the time of the solar UV-flux 
maximum, the 1 1 yr term is judged as in-phase with and pri- 
marily as a response to the solar UV-flux forcing. Regions 
where it is not in-phase indicate the presence of a decadal- 
scale, dynamical forcing or a nonperiodic, episodic forcing. 

The analyses follow a two-step approach to account for the 
effects of serial correlation at lag 1 in the data, following Tiao 
et al. (1990, Appendix A) and Remsberg (2008, Sect. 2.2). 
Initially, the MLR model terms are fit to the data time series 
for a given altitude and latitude bin. Next, the time series of 
the model/data residual is analyzed for its lag-1 autocorrela- 
tion coefficient, cp. Then, the model terms are reformulated 
to include the effects of <p, and the data time series are refit 
to get the final regression coefficients for each term. 

Several improvements for the SAGE II v7 data affect the 
results herein. The sunrise (SR) and sunset (SS) profiles are 
registered in altitude more accurately than before because 
of a correction in the use of the associated solar ephemeris 
data within the SAGE II algorithms (Damadeo et al., 2013, 
Fig. 2). Seasonal variations in the ozone are resolved more 
accurately, as a result. That revision has led to significant re- 
ductions in the MLR ozone residuals at low and middle lat- 
itudes. The (p values for the transformed terms in the MLR 
models are now larger for v7 than for v6.2, and they affect 
the coefficients of the analyzed trend terms and to a lesser 
extent the 1 1 yr terms. The limited duration, solar scan events 
of mid- 1993 through early 1994 and their associated SAGE 
II profiles are not part of the v7 data set. Some spuriously 
low, SR ozone data profiles are screened out for the low lat- 
itudes in September and November of 1991 because of ex- 
cessive extinction effects from the Pinatubo aerosol layer. 
A very few v7 ozone profiles are also not retained because 
their associated errors exceed 1 0 % within the altitude range 
of 50-27.5 km. Finally, continuous temperature time series 
from the Modern-Era Retrospective Analysis for Research 
and Applications (MERRA) system (Rienecker et al., 2008) 
are employed for the removal of Rayleigh scattering effects 
in the v7 algorithm. Rayleigh effects are a significant correc- 
tion for retrievals of the SAGE II ozone above about 55 km 
(Chuet al., 1989). 

The bin-averaged SR and SS data points have been assem- 
bled into a single time series, in order to decrease the time 
intervals between their observations at a latitude zone. This 
approach is helpful because having twice as many points that 
are separated by smaller intervals means that the semiannual 
and annual terms can be resolved better. However, the small, 
but systematic SR/SS differences of the ozone of the upper 
stratosphere and the mesosphere can lead to significant vari- 
ations for the de-seasonalized, time series residuals. In fact, 
those residuals are characterized by lag- 1 autocorrelation co- 
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Fig. 1. Time series of bin-averaged SAGE II SR and SS ozone num- 
ber density measurements (in cm -3 multiplied by 10 -11 ) at 25° S 
and 37.5 km (top) v6.2, (bottom) v7. Terms of the MLR model are 
indicated at the lower left, where SC refers to an 1 1 yr sinusoid term. 
The oscillating solid curve is the model fit to the data, while the 
straight line curve is the sum of just the constant and linear trend 
terms. 


efficients that are near zero or slightly negative. The approach 
used herein is to analyze the time series of SR points and then 
series of SS points. The constant terms of the MLR models 
are slightly different from the analyses of each series. Fol- 
lowing that, the SR points and the SS points of the separate 
time series are adjusted by half that difference before com- 
bining both sets of points into a single time series for the final 
MLR analysis. The lag-1 autocorrelation coefficients are less 
negative or more positive for the residuals of the final time 
series, as expected, since changes in zonal mean atmospheric 
ozone tend to occur slowly or over weeks rather than days. 

Figure 1 compares the v7 with the v6.2 ozone time se- 
ries for 25° S and 37.5 km that was analyzed previously by 
RL (their Fig. 2). The solid oscillatory curve is the fit to the 
points based on all the terms of the MLR model, while the 
straight line is the sum of just the constant and linear trend 
terms. In addition to the seasonal terms, one can clearly see 
that there is an 1 1 yr (or SC-like) term that is closely in-phase 
with the solar UV-flux maxima that occurred in 1991 and in 
2002. At this point it is noted that number density is plotted 
as increasing on the y axis in Fig. 1, even though it under- 
goes a nearly exponential decrease with altitude in the upper 
atmosphere. Nevertheless, when the UV flux reaches its max- 
imum in 1991 and 2002, Fig. 1 shows that ozone also attains 
its maximum (production) value at those times. A primary 
difference in Fig. 1 between the time series for v7 vs. v6.2 
is that the data points for v7 are defining the seasonal cycles 
more clearly; in general, the amplitudes of the annual and es- 
pecially the semiannual terms are larger with the analysis of 
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SAGE v7 03 1 1-yr Responses (%) for 91/05 



LATITUDE 


Fig. 2. Contour plot of the maximum minus minimum, 1 1 yr re- 
sponse (in percent) for the SAGE II v7 ozone of September 1991 
through August 2005. Contour interval is 1.0%. Shading denotes 
confidence intervals (Cl) > 90 % for the responses. 


the v7 data. That outcome is due to the more correct applica- 
tion of the solar ephemeris data for the registration of the SR 
and SS profiles. As a result, the residuals are smaller from 
the MLR fit to the v7 data, and the autocorrelation effects are 
larger. For instance, the coefficient <p is 0.44 for v7, but only 
0.05 for the v6.2 time series of Fig. 1. It was also realized 
that RL had inadvertently obtained bin-averaged points for 
their SAGE II ozone time series of 1991-2005 based on the 
poleward extremes for the latitudinal sweeps of the separate 
sunrise and sunset tangent track measurements, but using the 
dates for those extremes from the HALOE experiment rather 
than basing them on the actual beginning and ending dates 
for the sweeps of SAGE II. That incorrect grouping for the 
SAGE II profiles by RL also reduced slightly the accuracy of 
their fits for the seasonal terms. That error has been corrected 
for the results of Fig. 1 and for the remaining analyses herein, 
although it does not alter the conclusions of RL. 

Figures 2 and 3 represent an update of the MLR analy- 
sis results shown in RL for v6.2 (their Figs. 6-7), but based 
now on the v7 data. Figure 2 is the distribution of the 1 1 yr 
responses, and they are of the order of 3-5 % throughout 
most of the region with a minimum near 2 % at the tropi- 
cal stratopause. The analyzed responses increase to 6 % in 
the middle stratosphere of the northern subtropics, where Lee 
and Smith (2003) and RL indicated that there was an anoma- 
lous ozone forcing for some months following the Pinatubo 
event. A likely explanation for that ozone increase is from 
the effects of “de-noxification” or the temporary conversion 
of NO x to HNO3 in the presence of enhanced amounts of 
aerosols. Significance for the distribution of the amplitudes 
in Fig. 2 is indicated by the shading, which shows the do- 
main of the probability or the 90 % confidence interval (Cl) 


Phase of 1 1-yr max (in yrs) from Jan 1991 



LATITUDE 


Fig. 3. Contour plot of the phase variations (in years from January 
1991 or 2002) of the 1 1 yr response terms in Fig. 2. Contour interval 
is 1 yr. The phase domain of =b 1 yr is shaded and is considered as 
in-phase with the solar UV-flux maximum. 


that the 1 1 yr term is present in the zonally binned, ozone 
time series. Examples of further significance tests are given 
in Remsberg et al. (2001). Figure 3 is the companion plot of 
the phases of the 1 1 yr terms as referenced to January 1991, 
which is the approximate midpoint time of the maximum for 
the solar UV flux to the stratosphere. The dark shaded re- 
gions in Fig. 3 are where the phases are within ±lyr of 
January 1991 and where the 1 1 yr terms are considered as 
in-phase with the solar UV flux. Those 1 1 yr terms are inter- 
preted as the response to the UV-flux variations and are de- 
noted as SC-like. The 1 1 yr terms in the middle stratosphere 
of the northern subtropics are also closely in-phase in Fig. 3, 
even though they include perturbing atmospheric effects fol- 
lowing the Pinatubo event. 

Figure 4 shows the distribution of the coefficients of the 
associated linear trend terms for the SAGE II v7 ozone of 
1991-2005 in units of percent per decade, along with the 
shadings for their 90% and their 70-90% confidence in- 
tervals. Note that these trends are not explicitly related to 
an EESC term, but do include the linear effects of the rel- 
atively small changes in the atmospheric chlorine over this 
time span. There is much better continuity in the distribution 
of the trends in Fig. 4 compared to that from v6.2 in RL (not 
shown, but see their Fig. 8). This finding is expected, given 
that the accounting for autocorrelation effects in the residu- 
als has a primary effect on the trend terms. The trends from 
v7 are —4 % decade -1 from 40 to 50 km but —2 % decade -1 
from 30 to 35 km and also near 55 km. The anomalously 
large, negative trends above 50 km that were reported from 
the v6.2 data (see also RL, Fig. 8) are not present in the anal- 
ysis of the v7 ozone. This improvement may be because of 
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Linear Trends in SAGE 03 v7 from 91/05 (%/decade) 
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Fig. 4. Contour plot of the linear trend terms (% decade -1 ) from the 
MLR models for the SAGE II ozone data of 1991-2005. Contour 
interval is 2% decade -1 . Dashed contours denote negative trends, 
and the dotted contour is where the trends are zero. Darker shad- 
ing denotes where Cl >90% for the trends; lighter shading has 
70 %< Cl <90%. 



the use of the MERRA temperatures for the mesosphere in 
the v7 retrieval algorithm, particularly after 2000. 

Figure 5 displays the distribution of the SC-like response 
amplitudes for 1991-2005, according to the values that the 
1 1 yr terms would have on January 1991 . In other words, the 
1 1 yr amplitudes of Fig. 2 were adjusted for the fact that the 
phases of Fig. 3 were slightly lagging or leading the assigned 
time for the solar UV-flux maximum. Maximum responses 
of order 2. 5-4. 5 % occur in the upper stratosphere at low to 
middle latitudes, decreasing to 2-3 % at the higher latitudes. 
Although the SC-like responses are less than 2 % in the mid- 
dle stratosphere of the southern tropics, the phase-adjusted, 
SC-like responses of the northern subtropics are still anoma- 
lously large for this period because of the end-point anoma- 
lies from the Pinatubo event for the trend term of the MLR 
model fit. 

Figure 6 compares the SC-like ozone response results from 
SAGE II v7 and v6.2 for just the low latitudes. Specifically, 
the separate response profiles are averaged for each data set 
across the latitudes of 15° S, 5° S, 5° N, and 15° N. Because 
of the 20° width of each latitude bin, the average response 
profiles are characteristic of the latitude range of 25° S to 
25° N. There is little difference in the SC-like profiles from 
the two data versions. Figure 6 also shows the SC-like re- 
sponse profile obtained by RL for the HALOE ozone mix- 
ing ratio vs. pressure profiles, and they range from 0.5 to 
2.5 % for the period 1991-2005. The SC-like response from 
HALOE is smaller than that from SAGE II because ozone 
is more nearly conserved on pressure than on altitude sur- 
faces. The HALOE result indicates a small but increasing 
response in the mesosphere that is similar in shape to that 
from SAGE II, although neither is really significant. Fig- 
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Fig. 5. As in Fig. 2, but the SAGE II max minus min, 1 1 yr responses 
have been adjusted to be in-phase for that of a solar UV-flux maxi- 
mum in January 1991 . Shading denotes where Cl > 90 %. 


ure 6 also includes a low latitude response profile from the 
interactive, chemical-dynamical-radiative model of Brasseur 
(1993). There is good agreement between the response pro- 
files from HALOE and the model, both of which are given in 
terms of mixing ratio vs. pressure-altitude. They differ only 
in the lower mesosphere, where the response in the model is 
zero to slightly negative. 

Estimates of total 1 -a error are shown in Fig. 6 for the 
SAGE II v7 and the HALOE vl9 SC-like response profiles. 
The total error has two components. The first component is 
the uncertainty for the amplitude of the 1 1 yr response, which 
approximates very closely to the SC-like response in most 
instances. This component of the error increases with alti- 
tude in a manner that agrees with the doubling of the un- 
certainties from 3 to 6 % for the SAGE II ozone from 36 to 
39 km (Fioletov et al., 2006). Those errors are much larger 
than the amplitude of the SC response and reflect the effects 
of the decreasing signal to noise with altitude for the mea- 
sured SAGE II and HALOE profiles. In fact, the amplitudes 
of all the periodic terms are small at altitudes where ozone 
is in photochemical balance. Thus, the analyzed SC-like re- 
sponses from occultation data will not be significant for the 
uppermost stratosphere and into the mesosphere unless one 
analyzes data averages across many profiles. A second error 
component accounts for the variation of the response profiles 
from the separate latitude bins about the average of their re- 
sponses across 25° S-25° N. That component is large near 
30 km because of the perturbation in the retrieved ozone in 
the northern but not the southern subtropics following the 
Pinatubo event. This component of the error is larger for 
SAGE II than for HALOE because the measured signals for 
the SAGE II ozone are more sensitive to the interfering ex- 
tinction from the aerosols. However, both error components 
for SAGE II are essentially the same from the v7 and v6.2 
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SC-like 03 Response: SAGE, HALOE, Model 



Ozone Response. Max - Min (%) 


Fig. 6. Profiles of the average SC-like ozone responses (in per- 
cent) at the lower latitudes from the Brasseur model (dashed), the 
HALOE data (solid), and the SAGE II v6.2 (dotted) and v7 (dash- 
dot-dot-dot) data from 1991 to 2005. Horizontal bars are the RSS 
total errors for the HALOE and the SAGE II v7 ozone responses 
(see text). 

data (not plotted). Estimates of the total error are obtained by 
adding the two components using the root of sum of squares 
(RSS) method. It is also noted that the SC-like response pro- 
files of ozone in Fig. 6 are not as smooth for SAGE II vs. 
HALOE. Part of that difference is because the SAGE II mea- 
surements have a higher vertical resolution. In addition, at- 
mospheric ozone is not conserved as well at an altitude level 
(SAGE II) as on a pressure surface (HALOE). 


3 Effects of dynamical forcings and/or end-point 
anomalies 

Dynamical perturbations following the eruption of Pinatubo 
and/or end-point anomalies may be imparting biases for the 
analyses of linear trend terms from the time series of 1991 — 
2005 and to a lesser extent for their associated SC-like re- 
sponses. Two alternate time periods are evaluated for those 
effects: 1984-1998, which bridges the 1991 period; and late 
1992-2005, which eliminates an initial segment of the data 
following Pinatubo. 

3.1 Analyses for the period 1984-1998 

The SAGE II v7 results for 1984-1998 are shown in 
Figs. 7-9. Figure 7 is the distribution of the responses for the 
1 1 yr terms. The max minus min responses are on the order of 

4 % from 40 to 50 km but decrease to 2-3 % near 30 km, and 
there is no minimum response at the tropical stratopause like 
that of Fig. 2 for 1991-2005. An anomalous response of near 
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Fig. 7. As in Fig. 2, but the max minus min variations are from the 
SAGE II data from November 1984 through October 1998. 

5 % appears in a small region between 30 and 35 km at the 
Equator, but there is no indication of perturbing effects fol- 
lowing Pinatubo in the northern subtropics. Figure 8 shows 
the distribution of the phases of the 1 1 yr terms. Again, they 
are in-phase with the solar UV flux throughout most of the 
latitude/altitude domain, indicating that there are no signif- 
icant, out-of-phase dynamical forcings affecting the results 
in the upper stratosphere. However, at the Equator and be- 
tween 30 and 35 km the 1 1 yr term has its diagnosed maxi- 
mum closer to 1992-1993 than to January 1991. That phase 
lag of 1-2 yr exhibits good continuity across the latitude bins 
of the tropical middle stratosphere. The 1 1 yr terms are also 
not quite in-phase at the highest latitudes of Fig. 8, possibly 
due to the effects of decadal-scale, dynamical forcings. Fig- 
ure 9 shows the distribution of the coefficients of the linear 
trend terms, which are very similar to the results obtained 
by Wang et al. (2002) for the period of 1984-1999. The clear 
variation with latitude at 40-45 km is explained by the effects 
of changes in the reactive chlorine during this period (cf., 
Fleming et al., 2011, their Fig. 3). The anomalously large, 
negative trend of —8 % decade -1 and centered at 32.5 km and 
5° N in Fig. 9 is nearly offset in this 14 yr time span by the 
large positive, max minus min, response of 5 % (Fig. 7) that 
is not in-phase with the solar cycle (Fig. 8). Elsewhere, the 
succession of the trends from Fig. 9 and then Fig. 4 rein- 
forces the findings of others that the large declines in ozone 
of the 1980s were leveling out during the period 1991-2005 
and beyond (WMO, 2007, 2011; Kyrola, et al., 2013). 

Any perturbation to the ozone that occurs at the begin- 
ning of the time series of 1991-2005 will act as an end 
point anomaly for the determination of the linear trend term 
and will confound the coefficient of the 1 1 yr term. Fig- 
ure 7 shows that the region of rather large, 1 1 yr responses 
in the northern subtropics at 30-35 km are smaller than in 
Fig. 2, despite the fact that the data period following Pinatubo 
contributes to both time series distributions. As a further 
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LATITUDE 

Fig. 8. As in Fig. 3, but the phases of the 1 1 yr terms are from the 
SAGE II data from November 1984 through October 1998. 



LATITUDE 

Fig. 9. As in Fig. 4, but the trends are from the SAGE II data from 
November 1984 through October 1998. 

example, Fig. 10 is the ozone time series from 1984 to 1998 
at 5° N and 30 km. One can see that the MLR model terms fit 
the data well and that the Pinatubo event has less of an effect 
on the diagnosed linear trend term. 

The responses for 1984-1998 have a 1-2 yr phase lag in 
the tropical middle stratosphere, according to Fig. 8. It is in- 
ferred that some feature of the time series prior to 1991 is not 
being fit very well. An important consideration for the accep- 
tance of the terms of an MLR model is an examination of the 
time series of its residual to be sure that it contains no remain- 
ing periodic structure (e.g., Remsberg et al., 2001). As an ex- 
ample, the MLR fit to the ozone time series of 5° N and 30 km 
in Fig. 10 is systematically high in 1989-1990, followed by 
a slight low bias in 1991-1992. Those results suggest that 
the true amplitude of the QBO may have been larger than av- 
erage for that interval and that this term of the MLR model 
did not account for that interannual variation very well. As 
a result, the analyzed 1 1 yr term is influenced by the relative 
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Fig. 10. As in Fig. 1, but for the v7 ozone at 5° N and 30 km from 
SAGE II for 1984-1998. 


minimum in the observed ozone of 1989-1990, and the max- 
imum for that diagnosed term occurs in 1 992 rather than near 
the beginning of 1991. 

Hood et al. (2010) and Dhomse et al. (2011) tested for 
the likelihood of a separate perturbation from the El Nino- 
Southern Oscillation (ENSO) of 1988-1989 and its possible 
reinforcement of the amplitude of the QBO. Typically, ENSO 
anomalies at 16 km ascend to 30 km about 18 months later, 
or in this case, near to the time of the Pinatubo eruption in 
June 1991 (Garcia et al., 2007). Thus, Hood et al. (2010) and 
Dhomse et al. (2011) added a term to their regression mod- 
els related to a lagging ENSO index, but they found that it 
had very little effect at 30 km and above. Hood et al. (2010) 
and Chipperfield et al. (1994) also reported on a negative cor- 
relation between tropical ozone and NO 2 at 30-32 km from 
the effects of a vertical advection of NO y associated with the 
QBO forcing. Since Fig. 8 shows good coherence for the 1- 
2 yr phase lag of the maxima of the 1 1 yr terms across tropi- 
cal latitudes and from 30 to 35 km, it is presumed that those 
1 1 yr terms are confounded with the ENSO/QBO forcing and 
with the associated chemical/transport effects at that time. 

Figure 1 1 displays the distribution of the SC-like response 
amplitudes for 1984-1998, according to the values that the 
11 yr terms would be on January 1991. As for the period 
1991-2005 (Fig. 5), the llyr amplitudes of Fig. 7 are ad- 
justed for the fact that the phases in Fig. 8 are slightly lagging 
or leading the assigned time for the solar UV-flux maximum. 
Since the 11 yr response maximum is lagging January 1991 
by a year or more in the tropical middle stratosphere, its esti- 
mated SC-like response is reduced to less than 2 %. Brasseur 
(1993) and McCormack et al. (2007) obtained maximum re- 
sponses of only 2.5-3 % from about 35 to 40 km from the 
effects of the UV-flux variations. Analyzed responses of less 
than 2 % at the higher latitudes in Fig. 11 may also indicate 
effects of a chemical loss of ozone due to wintertime descent 
and mixing of NO x to lower latitudes that is enhanced at solar 
maximum (Garcia et al., 1984). 

Differences between Figs. 5 and 1 1 are found for the trop- 
ics between 45 and 50 km, where the analyzed SC-like re- 
sponses are larger for 1984-1998 than for late 1991-2005 by 
nearly a factor of 2. Solar flux values were also larger for the 
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Fig. 11. As in Fig. 5, but the SC-like responses are based on adjust- 
ments to Fig. 7 for 1984-1998. 

solar maximum centered in 1990-1991 compared with the 
maximum of 2002 (Smith and Matthes, 2008), but that ex- 
plains only part of the difference. Another part may be related 
to a bias in the analyzed ozone trends for one or both periods 
(cf., Figs. 4 and 9). Analyzed ozone responses are more eas- 
ily confounded with the trends for the ozone for the relatively 
short, 14yr time series of the present study. Since the upper 
atmospheric ozone decreases are larger in the 1980s than in 
the 1990s, it is likely that any confounding effects from the 
trend term are larger for 1984-1998. 

Were there any decadal- scale, dynamical forcings and as- 
sociated temperature variations that were in-phase with the 
solar UV-flux forcing during 1984-1998, or at least prior to 
September 1991? To check this possibility, time series of the 
SAGE II ozone at 5° S and 47.5 km are shown in Fig. 12 
for each of the two periods, 1984-1998 and 1991-2005. The 
data points reveal a pronounced amplitude minimum every 
year in January/February for the 1980s, but not from 1993 
through 1998 in the first series. Pronounced minima occur 
once again in 1999, 2001, 2003, and 2004 in the second se- 
ries. The first semiannual term in the MLR model provides a 
better fit to the time series of 1984-1998 than of 1991-2005. 
However, the analyzed amplitude of the semiannual term of 
the MLR model shows little to no change between the two 
time series; it is the amplitude of the annual term that is larger 
for 1984-1998 (4.3 %) than for 1991-2005 (3.4 %). As a re- 
sult, the MLR model residuals are larger for the 1991-2005 
series, and its associated SC-like terms are not being resolved 
as well. The irregular character of the annual ozone minima 
near the tropical stratopause indicates that they are episodic 
and short-lived. It is hypothesized that they occur due to sec- 
ondary, net circulations in response to extratropical gravity 
wave forcings (Shu et al., 2013) at times of polar winter 
warming activity (Hitchman and Leovy, 1986; Solomon et 
al., 1986). The ozone minima indicate a net downward trans- 
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Fig. 12. Time series of SAGE II v7 ozone at 5° S and 47.5 km for 
(top) 1984-1998 and (bottom) for 1991-2005. 


port of air having lower ozone values. Pawson and Naujokat 
(1999) reported that the Northern Hemisphere, winter polar 
vortex was unusually stable during successive years through 
the middle 1990s, which may account for the smaller annual 
ozone amplitudes during those years. Charlton and Polvani 
(2007) and Manney et al. (2005) also found evidence for 
more frequent, sudden stratospheric warming activity in the 
1980s and early in the 2000 decade, respectively. 

The average, SC-like response profile at low latitudes for 
the v7 ozone of 1984-1998 in terms of number density vs. 
altitude is compared in Fig. 13 with the profile for 1991— 
2005 that is replotted from Fig. 6. The ozone response from 
30 to 35 km is much smaller for 1984-1998 because there 
is less impact for the MLR analyses from the data just fol- 
lowing the Pinatubo event that occurs midway in this time 
series. The total i-o errors at 30 and 35 km are also reduced 
significantly and are more nearly equivalent to those from 
HALOE (cf., Fig. 6). However, the analyzed responses of 
Fig. 13 at 45 to 50 km do not differ significantly for the two 
SAGE II time periods. Dhomse et al. (201 1) conducted a re- 
gression analysis for the SC-like ozone response from their 
simulations of ozone based on the dynamical fields in the 
reanalysis products, ERA-40 and ERA-Interim (El), but for 
the longer time span of 1979-2005. They reported better dy- 
namical consistency in the El reanalysis, so their analyzed 
results based on El are considered here. Their analyzed SC- 
like response profile in terms of Dobson units (DU) per kilo- 
meter (or equivalent to number density) was shown in Fig. 6 
of their original discussion paper; that result is adopted and 
plotted in Fig. 13, too. Their simulated response profile has a 
peak value of 4 % near 40 km, and it agrees with the SAGE 
II v7 responses from 30 to 40 km. Presumably, the effects of 
any dynamical forcings due to an ENSO and/or following the 
Pinatubo event are present in the El product, and that is why 
the agreement is good from 30 to 40 km. Conversely, their 
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Fig. 13. As in Fig. 6, but the SC-like response profiles are for the 
model of Dhomse et al. (dashed; DU km -1 ) and from SAGE II v7 
for 1984-1998 (dash-dot) and for 1991-2005 (dash-dot-dot-dot). 
Horizontal bars are the RSS errors for the response profile from 
the SAGE II v7 ozone of 1984-1998. 


Fig. 14. As in Fig. 13, but the SC-like response profiles are for the 
model of Dhomse et al. (dashed in mixing ratio vs. altitude), from 
the HALOE (solid) and the SAGE II v7 (dash-dot) ozone for 1992- 
2005. Horizontal bars are the RSS errors for the response profile 
from the SAGE II v7 ozone of 1992-2005. 


analyzed SC responses decline sharply from 40 to 50 km and 
disagree with those from SAGE II. This altitude region is 
where both the trends in reactive chlorine and in temperature 
are important for characterizing the longer-term variations in 
ozone number density (see also Sect. 3.3). 

3.2 Analyses for the period late 1992-2005 

Lee and Smith (2003) showed that atmospheric forcing 
effects persisted for some months after the eruption of 
Pinatubo. Such forcings can lead to end-point anomalies for 
MLR analyses of ozone time series that begin by late 1991. 
The distribution of the SC-like amplitudes was obtained (not 
shown) for this period in the manner of the results in Figs. 5 
and 11. Then, Fig. 14 shows the analyzed SC-like response 
profile from SAGE II v7 for the low latitudes for the time 
series of late 1992-2005, or excluding the first 14 months 
of data following the Pinatubo event. An analogous response 
profile for 1992-2005 is shown for the HALOE ozone, but 
in terms of its mixing ratio vs. pressure values. The response 
from SAGE II varies from 3 to 4% between about 35 and 
40 km, decreasing to less than 2% at 30 km and at 50 km. 
The SAGE II response at 32.5 km for 1992-2005 (Fig. 14) 
is weaker, has a smaller uncertainty, and is more signifi- 
cant than for 1991-2005 (Fig. 6) because of the removal of 
the data points from the beginning of the time series. Both 
the SAGE II and HALOE responses appear to be free of 
the episodic forcing effects from ENSO and following from 
Pinatubo. The response at 42.5 km for 1992-2005 is also 
weaker than for 1991-2005, but the RSS errors show that 
it is not very significant. It may be that the rather large de- 


crease for the SAGE II responses at 45 to 50 km compared 
with those of 1984-1998 and 1991-2005 indicates that the 
11 yr terms of 1992-2005 are influenced by the lack of the 
Northern Hemisphere’s wintertime warming activity of the 
mid-1990s. 

An analyzed El-based ozone response profile for the 
1979-2005 time period that is in terms of mixing ratio (MR) 
vs. altitude is adopted from Dhomse et al. (2011, their Figs. 5 
and 6), and is included in Fig. 14. The El-based response 
profile includes the effects of the observed changes in the at- 
mospheric dynamics. Furthermore, the solar maximum, UV- 
flux values are larger for 1980 and 1991 than for 2002. Thus, 
their modeled response profile from the El-based time se- 
ries ought to be larger than that from the HALOE time series 
of 1992-2005. Figure 14 has a modeled response near 35 km 
with a peak value of 2.8 %, while the analyzed response from 
HALOE is 2%. Both response profiles decline to 1.5% at 
30 km and to near zero at 45 km, but then increase again from 
45 to 55 km. The shapes of the response profiles in Fig. 14 
agree well from SAGE II v7 and HALOE. It is concluded that 
the SC-like response profiles from SAGE II and HALOE are 
mainly due to variations in the UV-flux forcing for the period 
1992-2005. 

3.3 Effects of variations in temperature 

The SAGE II and HALOE responses in Fig. 14 differ in mag- 
nitude by 1-2 %, with the HALOE response being smaller. 
Vertical variations of atmospheric (and ozone) number den- 
sity are smaller on a pressure surface than at the correspond- 
ing altitude level (Akmaev et al., 2006). It is tempting to try 
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to relate the responses from SAGE II and HALOE by con- 
verting ozone number density vs. altitude (SAGE II units) 
to mixing ratio vs. pressure (HALOE units) or vice versa. 
However, the operational temperature data that accompany 
each of those satellite ozone data sets are uncertain (Randel 
et al., 2009). For instance, there is an apparent shift in the 
HALOE temperature time series in 2001 for altitudes below 
about 40 km (Remsberg and Deaver, 2005, their Fig. 4). This 
shift will lead to an incorrect trend of the pressure-altitudes 
that are used for the registration of the HALOE ozone pro- 
files. The MERRA temperature time series that are archived 
along with the SAGE II v7 ozone appear reasonable below 
40 km, but they display a separate shift in 1999 at 50 km and 
above. Nevertheless, it is emphasized that the analyzed ozone 
responses from SAGE II v7 are based on number density vs. 
altitude. Its algorithms do not depend on temperature, except 
for the removal of Rayleigh scattering effects that are signif- 
icant for the retrieved ozone in the middle mesosphere. 

Ozone is negatively correlated with temperature in the up- 
per stratosphere. Dhomse et al. (201 1) report a positive solar 
response in upper stratospheric temperatures, and it reduces 
their modeled SC ozone response signal. Are there periods 
when the observed, decadal-scale responses in temperature 
are trustworthy for the upper stratosphere? The analyzed SC- 
like temperature responses from two widely used reanalysis 
products, ERA-40 and ERA-Interim, were also provided by 
Dhomse et al. (2011); those max minus min responses var- 
ied from 1 to 2 K in the upper stratosphere. Remsberg (2009) 
also obtained an SC-like response for T (z) of about 1 K at the 
low latitudes at 40-45 km, which is an altitude region where 
temperature is a retrieved product from HALOE and where 
there are no apparent discontinuities in its time series. 

Radiative temperature changes due to the observed trends 
in ozone and CO 2 are largest at the low latitudes (Fleming 
et al., 2011). The rate of cooling in the upper stratosphere 
should have been greater in the 1980s, when ozone was be- 
ing depleted significantly (Akmaev et al., 2006). Were the 
observed trends in upper stratospheric temperature greater 
in the 1980s? Temperature time series data from several 
ground-based lidar stations are considered as a benchmark 
for changes in T (z), and they indicate that the cooling of the 
upper stratosphere was indeed larger for the period 1984- 
1998 than for 1992-2005 (Keckhut et al., 2011; Steinbrecht 
et al., 2009). Wang et al. (2012) and Thompson et al. (2012) 
also found global-mean, upper stratospheric temperature 
anomalies that are of the order of — 2 K decade -1 from 1984 
to 1995, but nearly constant thereafter. The HALOE experi- 
ment provides the only other near-global, satellite tempera- 
ture data set that extends for more than a decade during that 
latter period. Remsberg (2009) reported a trend in the upper- 
most stratosphere for HALOE T (z) that is of the order of 
only — 0.5 K decade -1 for 1991-2005. 

It is tentatively concluded that there is an SC-like response 
in T (z) of 0-1.5 K from 35 to 45 km, and that this response 
is very near to one from the El data used by Dhomse et 
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al. (2011, Fig. 5d) for simulations of their concurrent ozone 
response profile. It is also concluded that the trends in tem- 
perature, T (z), were significant during the 1980s, but nearly 
constant from 1992 to 2005. These changes in T(z ) along 
with the negative correlation between ozone and temperature 
explain qualitatively why there is a smaller SC-like response 
from the SAGE II ozone number density time series of 1 992- 
2005 (Fig. 14), as opposed to those of 1984-1998 (Fig. 13). It 
is also encouraging to see that the response profile in Fig. 13 
from the simulations of Dhomse et al. (201 1) for the longer 
period of 1979-2005 shows magnitudes that are midway be- 
tween those from the SAGE II v7 ozone from its two sepa- 
rate, shorter time periods. 

3.4 Ozone trends 

Ozone number density at an altitude is affected in the upper 
stratosphere by changes of the reactive chlorine. Figure 15 
shows the linear trends at 45 km from the analyses of the 
SAGE II v7 ozone for the two separate periods of 1984- 
1998 and 1992-2005. The vertical bars are the l -a uncer- 
tainties for the trend coefficients from the respective MLR 
analyses. The analyzed ozone trends for 1984-1998 for mid- 
dle latitudes of the Northern Hemisphere agree closely with 
the ozone decreases from —6 to — 9% decade -1 at 40 km 
reported by Steinbrecht et al. (2009) from ground-based li- 
dar measurements from 1979 to the late 1990s. They also 
agree closely with trends from the v7 ozone obtained by 
Kyrola et al. (2013). The trends from SAGE II are smaller 
for 1992-2005 or —3 to — 5% decade -1 at middle lati- 
tudes. At tropical latitudes the modeled ozone trends are 
from —3 to — 4% decade -1 and related to the changing at- 
mospheric chlorine (Fleming et al., 2011); the trends from 
SAGE II v7 support their model result. Tropical ozone trends 
from SAGE II v7 are less negative for 1992-2005 or about 
—0.5 to —2.5 % decade -1 and with 1 -a uncertainties of only 
zb 1 % decade - 1 . All other things being equal, this finding im- 
plies that reactive chlorine did not change by much between 
the two time periods. This finding agrees qualitatively with 
the average estimates of 2.8ppbv for total chlorine in the 
form of HC1 during the period 1984-1998 and of 3.1 ppbv 
for 1992-2005, or a difference of only 10 % (WMO, 2003). 

In general, the uncertainties for the trend terms are small 
and indicate that a linear term is a good fit to the data of 
the upper stratosphere for the relatively short, 13-14 yr time 
series shown in Fig. 15. Furthermore, even if the true trend 
for the time series is slightly nonlinear, that deficiency is ef- 
fectively accounted for by the inclusion of the lag-1 auto- 
correlation coefficient in the formulation of the terms of the 
final MLR models. Truncation of the ozone time series to 
minimize end-point anomalies following the Pinatubo event 
provides for significant trends from both the SAGE II and 
HALOE time series of 1992-2005. The small uncertainties 
for the results in Fig. 15 suggest that the analyzed trends 
are relatively unaffected by any bias in the associated 1 1 yr 
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Fig. 15. Ozone trends (% decade *) vs. latitude from SAGE II v7 
at 45 km from both 1984 to 1998 and late 1992 to 2005 and from 
HALOE at 1.7 hPa from 1992 to 2005. Vertical bars are the l -a 
errors for the coefficients of the trend terms. 


terms. Conversely though, a bias from the analyzed trend 
term can have a large effect on the diagnosed, SC-like ozone 
response profiles in Fig. 14. In fact, the amplitudes are small 
and sensitive to error for all the periodic terms of the MLR 
models in the upper stratosphere, as indicated for the SC- 
like response profiles in Figs. 6, 13, and 14. It is concluded 
that this circumstance can easily lead to disparities between 
analyzed SC-like response profiles from observed data time 
series vs. those obtained from model studies. 

The SAGE II results in Fig. 15 show negative ozone trends 
at 45 km that are larger in 1984-1998 than in 1992-2005. 
A part of that difference can be explained by cooling trends 
that are larger in the 1980s. In effect, ozone density should 
show a decrease with time due to a cooling or a contraction 
of its value at an altitude. Conversely, the HALOE trends in 
Fig. 15 are smaller than those from SAGE II for 1992-2005 
because the effect of the cooling is already accounted for 
in the analyzed trends for HALOE in terms of ozone mix- 
ing ratio on a pressure surface. There is also a hemispheric 
asymmetry in the ozone trends in Fig. 15 from SAGE II of 
up to 4 % decade -1 for 1992-2005, a result that is in accord 
with findings of Kyrola et al. (2013). The ozone trends from 
HALOE at 1 .7 hPa are more nearly symmetric across the two 
hemispheres. The asymmetry in the SAGE II ozone is most 
likely a result of dynamical forcings. For instance, Fig. 16 
shows the trend curves with latitude at the constant altitude 
of 45 km from SAGE II along with the concurrent SC-like 
response variations for the two time periods. There is hemi- 
spheric symmetry for both the SC-like response and the trend 
curves for 1984-1998, as expected from zonally averaged 
photochemical models. However, one can see that for 1992- 
2005 the SC-like responses increase positively with latitude 
in the Northern Hemisphere, while the trends increase neg- 
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Fig. 16. Trends and SC-like responses vs. latitude for SAGE II 
ozone at 45 km for the periods of 1984-1998 and 1992-2005. Ver- 
tical bars on the trend curves are their 1-cr errors. 

atively. A similar result is obtained from analyses of the 
HALOE ozone at the constant altitude of 45 km (not shown). 
Midwinter, wave-induced warmings were significant for the 
Northern Hemisphere for those years near solar maximum 
and lacking at solar minimum during the period 1992-2005, 
leading to a decadal-scale, vertical variation of the ozone sur- 
faces near 45 km that was enhanced and in-phase with the so- 
lar forcing (e.g., Shu et al., 2013; Gerber, 2012; Solomon et 
al., 1986; Pawson andNaujokat, 1999). A further assessment 
of this likelihood is beyond the scope of the current study. 
Nevertheless, it is judged that the differences for the SAGE II 
and the HALOE ozone trends in Fig. 15 and for the analyzed 
SC-like responses for SAGE II in Fig. 16 are trustworthy for 
1992-2005 because atmospheric perturbations that followed 
from the Pinatubo event have very little effect on the ozone 
time series of 1992-2005. 


4 Comparisons with other published SC-like ozone 
response profiles 

SC response profiles have been reported by RW (their 
Fig. 12a) based on their analyses of the SAGE I plus SAGE 
II v6.2 or the so-called SAGE-based ozone. They analyzed 
ozone time series at altitudes in terms of integrated ozone or 
DU km -1 , which is equivalent to ozone number density vs. 
altitude. Their response profiles for 1979-2005 have maxi- 
mum values of about 4 % at middle latitudes (after adjust- 
ing their results by the actual units of the FI 0.7 flux maxi- 
mum), and the distributions herein (see Figs. 5 and 11) agree 
with them. However, their responses at the low latitudes are 
smaller, most likely because they are unable to fully resolve 
the response maxima from v6.2 and their relatively large in- 
terannual residuals compared with those from v7. Their in- 
clusion of the SAGE I ozone may also be introducing biases 
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for their analyzed ozone trend terms; the SAGE I ozone has 
yet to be reprocessed through the current v7 algorithm. RW 
considered monthly averages of SAGE profiles at every 5° of 
latitude, or at much narrower bin widths than for the present 
study. Thus, it is expected that the random uncertainties for 
their analyzed SC-like terms ought to be larger by a factor of 
two. In other respects the component of the v7 ozone error in 
Figs. 6, 13, and 14 that increases with altitude is qualitatively 
similar to the levels of unexplained, interannual variance re- 
ported by RW (their Fig. 3). 

Hassler et al. (2008) updated the ozone time series and 
climatology of RW, and merged them with ozone data from 
other, more recent satellite and ozonesonde measurements. 
McLinden et al. (2009) assembled a SAGE-corrected, time 
series of SBUV ozone for the purpose of providing better 
estimates of its long-term trends. Dhomse et al. (2011) an- 
alyzed those same SAGE-corrected, SBUV ozone data and 
compared those results with the findings from their simula- 
tion studies. Still, all of their studies used the SAGE v6.2 
data, and the adjustments to the SBUV data of 1978-1981 
were based on the SAGE I ozone that has yet to be repro- 
cessed. 

SH also analyzed SAGE II v6.2, HALOE vl9, and SBUV 
v8 ozone for the time span of 1979-2003, but after de- 
seasonalizing the data by removing their respective clima- 
tological mean values. For their analyses of the time series 
of the 3 -monthly averages of the SAGE data, they did not in- 
clude ozone from SAGE I nor from SAGE II past 2003; they 
also excluded just over two years of data following the erup- 
tion of Pinatubo. The present analyzed responses are similar 
to theirs for the middle latitudes. However, SH obtained max 
minus min responses for the low latitudes that are less than 
2 % between 28 and 36 km. The present study finds responses 
from SAGE II that are also small at 30 km, but which then in- 
crease to 3 % at 35 km (see Figs. 13 and 14). It is judged that 
at least part of the difference with SH is a result of the im- 
proved accounting for the effects of the seasonal cycles in the 
SAGE II v7 data. 

SH and Toupali et al. (2007) analyzed time series of 
monthly SBUV v8 ozone profiles from 1979 to 2003 for their 
solar cycle responses. Their analyses for the SC responses 
were performed on separate time series for each of the in- 
dividual months or seasons. Both studies reported near zero 
SC responses from about 30 to 40 km at the low latitudes, 
contrary to the predictions of most photochemical models. 
Their response profiles from the SBUV ozone mimic those 
shown by Torres and Bhartia (1995, their Fig. 4) for the re- 
trieval effects on ozone at higher altitudes, due to not ac- 
counting for multiply- scattered, UV channel radiances from 
a volcanic aerosol layer located at 26 km. The effect of the 
aerosol scattering in the retrieval in that instance is a signif- 
icant reduction in the ozone from 26 to 40 km and a small 
increase above 40 km. Since the Pinatubo aerosol layer was 
present near solar-UV maximum, it is very likely that the ef- 
fects of the scattering completely canceled the effects of the 
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enhanced solar-UV flux. It is noted that SH did not exclude 
any SBUV data following the volcanic events of El Chichon 
and Pinatubo for their analyses, while Tourpali et al. (2007) 
deleted data for 1982 and 1991. RW were more conservative 
for their regression studies; they deleted two years of their 
merged SBUV v8 ozone following both volcanic events. 

The SC-like responses of SH from the SBUV and SAGE 
II data are included in WMO (2007, Figs. 3-19) along with 
simulated responses from several models. One minor puzzle 
has been the increasing SC response values from the obser- 
vations of SAGE II v7 and of HALOE at 50 km and above 
(e.g., Figs. 13 and 14), even though they are not very signifi- 
cant from the present analyses. Responses from the Brasseur 
model (and from many others) are near zero to negative at 
these altitudes. Yet, the response from the 3-D model study 
of Tourpali et al. (2003) is of the order of 2% at 50 km, 
and that from the simulation study of Dhomse et al. (201 1) 
increases from near zero at 45 km to 2% at about 55 km. 
It may be that the increasing responses with altitude above 
the stratopause from the SAGE II and SBUV ozone anal- 
yses are largely real and a reflection of dynamical forcing 
mechanisms. As an example, Tourpali et al. (2003, Fig. 1) 
obtained SC-like ozone responses of 2 % or more throughout 
the winter hemisphere at and above 55 km, presumably due 
to dynamical rather than direct solar-UV forcings. A mid- 
to-late winter stratospheric warming and its associated lower 
mesospheric cooling activity can modify the annual ozone 
amplitude. Of course, such events tend to be more episodic, 
such that an MLR analysis would not fit the annual varia- 
tions in the observed ozone variations well every winter. If 
mid- winter warming activity occurs more frequently near so- 
lar maximum than solar minimum (or vice versa), then it can 
represent an additional, decadal-scale forcing for the ozone 
in the upper stratosphere/lower mesosphere. 

The role of changes in water vapor on ozone should not 
be ruled out either, especially for the mesosphere. Relatively 
low values of H 2 O will mean less ozone destruction via the 
HO x chemical cycle. Nedoluha et al. (2009) showed from the 
HALOE data that water vapor in the lower mesosphere was 
increasing at a rate of 1-2 %yr -1 from 1991 through 1996 
and then remained fairly steady to 2005. Such a trend will 
affect the analyzed solar response and trends in the ozone for 
1992-2005. Less H 2 O (and HO x ) near solar max in 1991- 
1992 will mean that there is more ozone at chemical equi- 
librium with all other factors being constant. In situ observa- 
tions in the lower stratosphere indicate that H 2 O was also in- 
creasing at a fairly steady rate of 1 % yr -1 through the 1980s 
(Rosenlof et al., 2001). It is expected that those positive 
trends also extended to the mesosphere during that decade, 
but with a lag of several years. Thus, for the relatively short 
ozone time series of the present study the effects of such 
trends in H 2 O can confound an analysis of the concurrent 
SC-like term. More detailed studies of the decadal-scale ef- 
fects of mid-winter warming events and/or trends in water 
vapor are beyond the scope of this report. 
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5 Conclusions 

This study is an extension of the analyses and findings in 
RL, but it is based now on the newly processed SAGE II v7 
ozone. MLR analyses for the period 1991-2005 were con- 
ducted, as before, and then additional analyses were made 
for the SAGE II v7 ozone of 1984-1998 and 1992-2005. 
The seasonal cycles fit the ozone time series of v7 much bet- 
ter than for the previous v6.2, and the MLR de-seasonalized 
residuals are smaller throughout much of the stratosphere for 
v7. As a result, the interannual and trend terms are resolved 
better and have greater significance than before, especially 
from 27.5 to 42.5 km. It is judged that the SAGE II v7 ozone 
is of better quality. 

The max minus min, 1 1 yr ozone responses from SAGE II 
v7 for the period of 1984-1998 differ somewhat from those 
of 1991-2005 in the middle and upper stratosphere, but they 
are in-phase with the solar cycle UV flux at most altitudes 
and latitudes. The shapes of the SC-like response profiles 
from the SAGE II and HALOE ozone agree for the period 
of 1992-2005, or after deleting the first 14 months of data 
following the Pinatubo event and after taking into account 
the geophysical units of their respective data sets. At tropical 
latitudes and between 30 and 35 km the 1 1 yr responses from 
SAGE II for 1984-1998 were lagging the UV-flux maximum 
by 1-2 yr because of a separate perturbation that affected its 
ozone time series during the period 1989-1990. 

Ozone trends in the upper stratosphere from SAGE II v7 
extend to — 8% decade -1 at the middle latitudes for 1984- 
1998. Model studies of the effects of atmospheric chlorine 
are generally consistent with that finding. The trends for 
1992-2005 vary from about —4% decade -1 to near zero 
from the Northern to the Southern Hemisphere. Differences 
between the ozone trends from SAGE II at 45 km and from 
HALOE at 1.7 hPa are also about 4% decade -1 at northern 
middle latitudes but near zero at southern latitudes. Based 
on published diagnostic studies of others, it is inferred that 
the effects of wintertime wave forcings and their associated, 
induced net circulations are contributing to the difference in 
the Northern Hemisphere. 
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